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The DAMA dark matter search experiment observes a statistically significant percent-level 
variation of its low-energy count rate with a period of one year. In this note we recall some of 
the arguments which challenge the hypothesis that the cosmic ray induced underground muon 
flux can be the origin of the modulation. In addition, we provide new comments on recent 
■ works on this subject. 

o 

^ ' 1 Introduction 

>v 

In a recent publication^ henceforth referred to as CPY we examine recent claims that the cos mic 
ray muon flux can be responsible for generating the modulation signals seen by the 

damaEE] 

and, more recently, by the CoGeNT™ Dark Matter (DM) direct detection experiments. By 
examining the time series of the reported DAMA and CoGeNT signals we find the data sets 
differ in phase, likely in amplitude, and potentially in their power spectrum from measurements 
i-C . of the underground muon flux. A correlation analysis in CPY reveals that the data sets are 

incompatible and that the muon flux as the sole source of the reported signals is excluded at 
^■r high confidence. 

In this note we restrict the discussion of the muon hypothesis to DAMA which is the only 
direct detection experiment which claims detection of a firm DM signal. It is situated at Gran 
Sasso, Italy, in the underground LNGS laboratory. The collected data encompasses two major 
runs, DAMA/Nal (Dec 1995 - July 2002) and DAMA/LIBRA (Sept 2003 - Sept 2009). The 
residual count rate reported by the collaboration exhibits annual modulation compatible with 
\Q ' what is expected from generic DM models. 

TeV-scale cosmic ray muons produced at stratospheric altitude levels can reach deep under- 
ground and induce spallation reactions in the detector and nearby. The seasonal variations of 
the underground muon flux have been observed on the northern and southern hemisphere alike. 
For DAMA, the relevant measurement at the LNGS site is the one from LVD^l which over- 
laps in time with DAMA/LIBRA cycles 1-5. The integral muon intensity underground reads 
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(In) — 3 x 10~ 4 mT 2 sec -1 and exhibits an annual variation of ~ 2% in amplitude. In Fig. Q] 
we show the percent residuals of the muon flux when binned in concordance with DAMA with 
the annual mean count rate subtracted. In addition, the (2-4) keVee bin of the DAMA/LIBRA 
experiment are shown; a baseline rate of s = 1.15 cpd/kg/keV has been assumed. The seeming 
similarity in time and amplitude is tantalizing and has lead to the n otion t hat both signals may 
in fact be measurements of the very same cosmic ray phenomenon EH3 i n the following we 
recap and further develop on the arguments which challenge this hypothesis. 

2 The Null Hypothesis: periods, binning, and subtraction of means 

A spectral analysis will most readily reveal periodic behavior in a time series {tj, ^Q}»=i at- br 
case the times ti are unevenly spaced, it is most convenient to make use of the Lomb-Scargle 
(LS) periodogram^, 
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Here r is obtained from tan(2o;r) = ]T\ sin2o;ij/E i cos 2wij, and it ensures that Px{oS) is 
invariant to shifts in the time origin; X and a\ are the arithmetic mean and the variance of 
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Figure 1: Percent annual residuals of the LVD measured muon flux when binned in accordance with the 
DAMA/LIBRA runs 1-5. The latter residuals are shown for the 2-4 keV bins assuming a baseline count rate of 

s= 1.15 cpd/kg/keV. 



the data {Aj}, respectively. A Bayesian derivation and generalization of (pQ) can be found in 
CPY. Noise is rejected at the 1 — a confidence level by demanding that the power at a sampled 
frequency is in excess of 
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A statistical penalty has been included for examining M independent frequencies. In order to 
obtain a meaningful statistical interpretation we have to choose M carefully. In CPY we chose to 
sample the data at frequencies w n = nujp with n = 1, . . . , N, M = N, and with the fundamental 
frequency ujp = 2tt/T where T is roughly the range of years covered by the data set. In this 
note we instead oversample the data using a finer grid of frequencies. Since Xj are not grossly 
unevenly spaced in tj, M ~ N holds^l a g long as we only sample Aj above the average Nyquist 
frequency, v c = 1/2 A where A denotes the average spacing between tj. We checked that both 
approaches agree on coinciding frequencies. 

The left panel of Fig. [2]shows the DAMA/LIBRA power spectrum for the residual count rates 
as reported in^ and which are obtained by subtracting the mean count rate in each cycle^^ for 
the various energy bins as labeled. Above the horizontal line, noise is rejected at 99% C.L. With 
an approximate power Px(27r/lyr) ~ 12 at a period of 1 yr, the null hypothesis is disfavored 
by at least 4.5<r at a single frequency and at 3.6<7 according to ([2|). Notably, there is significant 
power at a period at 480 days in the (2 — 4) keV ee and (2 — 6) keV ee bins and a consistent power 
in all three energy bins with period 1 /3 yr corresponding to a triannual mode; for the latter only 
the (2 - 5) keV ee bin exceeds the 99% C.L. level. 

The right panel of Fig. [2] shows the LS periodogram obtained from LVD data where only 
those data points were included which have actual temporal overlap with the bins of the 
DAMA/LIBRA cycles, namely, full cycles 1-5. We only show a fraction of the full power 
spectrum with smaller periods omitted. The horizontal line is calculated according to ([2]) from 
all sampled independent frequencies. Since Alvd 3> Adama ^ e rec l u i r ed power to reject the 
null hypothesis is larger when compared to the right plot — despite the similar period intervals 
shown. The solid line is obtained from the unbinned LVD data set. Aside from the significant 
yearly modulation, the spectrum exhibits temporal variations on time-scales smaller and greater 
than one year. At first sight this is in stark contrast to the DAMA spectra which barely exhibit 
power at periods other than at T = 1 yr. However, in order to better compare these data sets 
one should (1) bin the muon data over the same time spans DAMA/LIBRA bins their signal 
and (2) subtract the yearly mean of muon flux. The latter mimics the procedure employed by 
the DAMA collaboration when they present their residual count rates. To isolate the effects of 
(1) and (2) we carry out each procedure separately and show its effect on the power spectrum. 




Figure 2: Left: LS spectrum for DAMA/LIBRA residual count rates as a function of the period. The horizontal 
gray line marks the power above which noise can be rejected with 99% confidence for any frequency. Right: 
Power spectra of the LVD data which has temporal overlap with DAMA/LIBRA (cycles 1-5). The solid line is 
for unbinned LVD data; the dotted line is obtained after binning the data in concordance with the time bins of 
DAMA/LIBRA and after re-scaling the resulting power spectrum to match significance levels of the solid curve: 
the dotted line is no LS but can directly be compared with the other lines. The dashed line is obtained from 
unbinned LVD data but the yearly mean intensity has been subtracted. The subtraction results in damping of 
the long period modes; binning damps power on all time-scales. Both effects make the muon data appear more 

similar to DAM A. 



We first observe that binning reduces the number of data points from N ~ 1.3k to N = 38 
when considering DAMA/LIBRA cycles 1-5. For a given significance a this changes the condi- 
tion on the minimum power in Eq. ([2]). Therefore, we first bin the LVD data and subsequently 
scale the resulting power spectrum in such a way that for a given y- value the significance agrees 
with the one from the unbinned data. The result is shown by the dotted line in Fig. [2j Though 
the dotted line is not a LS anymore, one can now directly compare the dotted with solid line. 
We observe that the generic effect of binning is to mask power on all time scales. The effect 
of subtracting the mean intensity from the data on a yearly basis is shown by the dashed line 
instead. This time the power spectrum is not scaled and the dashed line is indeed a LS peri- 
odogram. In contrast to binning, this procedure largely preserves power on periods T < lyr. 
However, for larger periods we again observe that variations are masked due to the particular 
treatment of the data. 

In contrast to the LVD spectrum, the Fourier transformed DAMA residuals in the left panel 
of Fig. [2] show much less significant power at periods smaller and greater than a year. As 
shown above, this may simply be an artifact of the way the DAMA/LIBRA modulation signal 
is obtained. It is unfortunate since a Fourier analysis could otherwise be used to discriminate 
between background induced effects and real signal. It is not clear how to assign a quantitative 
significance to this difference without knowledge of the full, unsubtracted time series. This calls 
for a release of the yearly baseline count rates by the DAMA collaboration. 

The DAMA collaboration has recently responded to the criticism that baselines are sub- 
tracted per cycle^21. Instead of providing the actual values, the LS power spectrum of the joint 
DAMA/Nal and DAMA/LIBRA baselines is presented. No significant power which would point 
towards long-term behavior is observed. However, the baselines are a mere set of 13 numbers. 



It is not surprising that no statistically significant power is observed from such a small data 
set. Similarly, we find that the LVD baseline spectrum does not yield any significant power 
when calculated according to (H|), too. This casts doubt that a presentation of the baseline-LS 
precludes the existence of variation s in D AMA with periods larger than one year, similarly to 
what is observed in muons; see also™. 

3 The phase of muons: a potentially stretched concept 

Given an isotropic DM velocity distribution f(v) in the earth's vicinity, the differential recoil 
rate cLR/cLEr of DM scattering in the detector is predicted to peak on June 2 nd , corresponding 
to a phase of to = 152.5 days after January 1 st , 

dR [°° f(v) , r , 

oc / dv Cq + c\ cos [u[t — to)\. (3) 
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Here, v m i n is the minimum required relative velocity between DM and target which produces 
a nuclear recoil of energy Er. In the canonical case, the annual mode is at the percent level, 
ci/co = C(10 -2 ); higher harmonic corrections to ([3]) have first been pointed out in CPY. The 
clear-cut prediction on to arises from the geometric setup of the earth's rotation around the sun 
in conjunction with the movement of the solar system relative to the Galactic coordinate sys- 
tem. Hence, apart from the period, the second most important characteristic of the oscillations 
observed by the DAMA experiment is the phase of the signal. 

On the other hand, when considering environmental factors as potential background, the 
variation in time will in general not be modulated with a strict period of a year. Additionally, 
we cannot expect a sinusoidal behavior similar to ([3]). Both assertions can be tested on the 
LVD and DAMA data sets. For this we model the data by a sinusoid, minimize the usual x 2 
function, and assess the goodness-of-fit. Confidence regions in to and T can be inferred from a 
Ax 2 -method which is equivalent to constructing a profile likelihood. Fixing the period to one 
year we find that the DAMA/LIBRA (2 - 4) keV ee and the full LVD data do not agree with 
respective values 

DAMA/LIBRA : to = (130 ± 8) days, x 2 /dof = 37.5/41, (4) 
full LVD : t = (185 ± 1.5) days, x 2 /dof = 3450/1971. (5) 

The uncertainty quoted is statistical. An additional uncertainty of ±2 days may be attributed 
to digitization. A few comments are in order: 1) The sinusoid provides a very poor description 
of the LVD data with prohibitively large x 2 ■ This supports our expectation that a complex 
phenomenon such as the underground muon flux cannot be adequately described by this simple 
model. 2) Even if we scale up the error bars for LVD to the point where x 2 /dof ~ 1, the 
stringent error on the phase remains below ±2 days. 3) An error of ±15 days on to as quoted 
by LVD0 is certainly incorrect. 4) Once T is allowed to float, the confidence regions in to and T 
become sensitive to the time-ori gin: for further details we refer the reader to CPY. These results 
have recently been confirmed in^^. One major shortcoming of this approach is that the data 
has to be subjected to the simple model of sinusoidal variation, before any conclusions can be 
drawn. Indeed, the power-spectrum of the LVD data makes it clear that significant power exists 
in modes with periods larger and smaller than one year. The notion of a single phase must 
therefore be treated with caution. 

4 A correlation analysis: the conclusive approach 



A glance at Fig. [T] suggests that the DAMA/LIBRA and LVD data sets are correlated. Indeed, 
when evaluating Pearson's coefficient of correlation, r € [—1,1], one finds for the DAMA and 



LVD sets as presented in Fig. Q] 



r LVD,DAMA = 0.44. (6) 

This value excludes the no-correlation hypothesis with a confidence level greater than 99%. How- 
ever, just because the data sets are correlated does not imply that they are causally connected. 
This can only be answered on a model- by-model basis. 

Clearly, a muon background model for DAMA must encompass the stochastic nature of 
the underlying process. Even if the timing between muons and DAMA are incommensurate 

at first sight, could it be that such a Poisson smearing alleviates the observed tension in the 
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annual phase? That this is indeed the case has been shown in 121 where a simple, linear model for 
how a muon-sourced background may be realized is employed. Taking from this that the muon 
hypothesis for DAMA would become viable again is however not correct. 

In CPY we generate 10 4 realizations for DAMA based on the generic model presented in^ 
which connects the muon flux measured in LVD with the count rate observed in DAMA/LIBRA. 
The question which has to be asked is now the following: How likely is it that one realization 
out of the mock data indeed induces the signal observed in DAMA. To answer this question 
we can again look at the correlation coefficient. It has the merit that we do not impose any 
functional form on the data. As expected, the generated mock data exhibits a high degree of 
correlation with LVD. In fact, the correlation with the actual muon data is substantially higher 
than the one which is actually observed in ([6]). This rules out the model at 99% C.L.. Given the 
generality of the model, it also strongly disfavors the hypothesis that muons can be the origin 
for the count rate variation seen in the DAMA detectors. 

5 Conclusions 

In CPY we lay out a detailed time-series analysis of DM direct detection data as well as related 
datasets. In particular, we find no evidence to support the claim that atmospheric muons are 
responsible for the signal that the DAMA collaboration observes. We identify difficulties in 
phase, amplitude, power spectrum and degree of implied correlation between the data sets. 
Here we offer some additional comments: 

• When binning data, fine-grained information on the time-series is lost. In addition, we find 
that when the LVD measured muon flux is binned in time in the same way as DAMA/LIBRA 
presents their data, significant power is lost on all time scales. This has the consequence that 
the LVD and DAMA data sets become more alike in the frequency domain with a dominating 
peak at a period of one year. Thus binning seems to diminish discriminating potential between 
the data sets. This calls for an unbinned analysis of the DAMA data set. 

• The uncertainty on the phase of the underground muon flux reported by LVD, to = (185 ± 15) 
days is erroneous. In CPY, an uncertainty of ±2 d ays was obtained, but the discrepancy had 
not been emphasized. An independent analysis iiv^l confirms our value. 

• A recent paper by the DAMA collaboration presents the LS spectra of the average count rates 
in the DAMA/Nal and DAMA/LIBRA cycles. No significant power indicating potential long- 
term variation (T > 1 yr) is observed. Following this procedure for LVD data, we can reach 
analogous conclusions — in seeming contradiction with the power observed in the right panel 
of Fig. [21 One cannot expect statistical significance from a LS-transform of O(10) numbers 
without further insight into the normalization of (HJ. 

• The same paper confirms our previous findings in CPY that a putative cosmic ray activation 
can only lead to a maximum phase shift of a quarter year in the subsequently modulated 



decay. Hence it is not possible to overcome the ~ 11 month discrepancy between muons and 
DAM A. 



• A recent paper^ incorporates additional data on the Gran Sasso underground muon flux. 
Prequentist fits to period and phase improve on the already significantly disjoint LVD and 
DAMA regions presented in CPY. This corroborates the statement that both data sets are 
seemingly incompatible. However, as argued in Sec. and at great length in CPY reliance 
on the inferred phase can be misleading. Indeed in^ and CPY it was shown that Poisson 
smearing can sufficiently alleviate the tension between the phases of the muon and DAMA 
data sets. It is the correlation analysis which reveals that the data sets remain nevertheless 
incompatible. However, the latter conclusion depends on the model which connects the muon 
flux with the DAMA count rate and it can only be achieved considering data sets which have 
actual temporal overlap. From this perspective, no further gain of information with respect 
to DAMA is obtained by enlarging the set of muon measurements as done in^l 

As in the case of DAMA, in CPY — employing a similar line of analyses — we find no significant 
correlation between the CoGeNT data set and the yearly variation of the muon flux. In sum- 
mary, cosmic ray muons cannot be the sole source for the observed modulation signals in these 
experiments. 
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